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Abstract 

We present an improved measurement of the rj' meson energy spectrum in T(1S) decays, using 
1.2 fb _1 of data taken at the T(15) center-of-mass energy with the CLEO III detector. We 
compare our results with models of the rj gluonic form factor that have been suggested to explain 
the unexpectedly large B — > n'X s rate. Models based on perturbative QCD fail to fit the data for 
large rf energies, and thus an explanation outside the realm of the Standard Model or an improved 
understanding of non-perturbative QCD effects may be needed to account for this large rate. 

PACS numbers: 13.25.Gv,13.25.Hw,13.66.Bc 
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I. INTRODUCTION 



CLEO observed a surprisingly large branching fraction for the decay B — > rfX s at large 
momenta of the rf meson p v > : BiB — > r/Xg)^ />2 Gev = 

BaBar ^ later obtained B(B -> r/'X s )| Pv >2 GeV = (3.9 ± 0.8 ± 0.5 ± 0.8) x 10" 4 . This rf 
momentum is beyond the end point for most b — > c decays, so the rf yield from 6 — > c 



is expected to be only of the order of 1 x 10 4 . Predictions assuming factorization 



lIU 



estimate the charmless component of this branching fraction to be also about 1 x 10 4 . 
Thus conventional calculations cannot account for the measured rate and they also fail 



to predict the right shape for the rf momentum spectrum [6(. These findings motivated 
intense theoretical activity because new physics could account for such an enhancement. 
However, Standard Model explanations have also been proposed. For example, Atwood and 
Soni [7j argued that the observed excess is due to an enhanced b — > sg penguin diagram, 
complemented by a strong rfgg* coupling, induced by the gluonic content of the rf wave 
function. Fig. [T] (left) shows the corresponding Feynman diagram. The high q 2 region of the 
(fgrf vertex function involved in this process also affects fast rf production in T(IS') decay 
[6, 7, sj], whose relevant diagram is shown in Fig. [T] (right). Thus a precise measurement of the 
rf inclusive spectra from the process T(IS') — > ggg* — > rfX can improve our understanding 
of important B meson decays. 
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grf) (left) and T(15) -» ggg" -> rfX (right). 



FIG. 1: Feynman Diagram for b — > s{ 



The effective vertex function rfgg* can be written as 7| H(q 2 )S ab e a /3 fMl/ q a k l3 eie2, where q is 
the (g*) virtual gluon's four-momentum, k is the (g) "on-shell" gluon's momentum (k 2 = 0), 
a, b represent color indices, £x; e 2 are ^ ne polarization vectors of the two gluons, and H(q 2 ) 
is the g*grf transition form factor. Different assumptions on the form factor dependence 
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have been proposed 



While ARGUS was the first experiment to study the inclusive rf production at the T(15') 



12j, they did not have enough data to separate T(IS') — > ggg* from the other components 



discussed below. CLEO II 13j was the first experiment to have sufficient statistics to measure 
inclusive 7/ production from the subprocess T(IS') — > ggg*. These data ruled out a class 



of form factors characterized 
theoretical calculations 



)y a very weak q 2 dependence Subsequently, several 

ll| derived the perturbative QCD form factors from models 



of the 7/ wave function. Attempts to use CLEO II data to constrain the model parameters 

n 

|14j were not conclusive, due to the limited statistics at the end point of the 77 spectrum. 
Thus, it was difficult to establish whether neglecting higher order terms in the perturbative 
expansion was appropriate ljj]. An improved measurement, based on a higher-statistics 
sample, is important to provide an improved determination of the QCD parameters, and, 
consequently, a more stringent test of the theory. This work reports a new measurement 
of the inclusive 7/ spectrum from the process T(1S) — > ggg* — > rfX based on the largest 
Tflff) sample presently available, more than a factor of 11 greater than the previous study 
13]. 



II. DATA SAMPLE AND ANALYSIS METHOD 

We use 1.2 fb _1 of CLEO III data recorded at the T(IS') resonance, at 9.46 GeV center- 
of-mass energy, containing 21.2 x 10 6 events, and off-resonance continuum data collected at 
center-of-mass energies of 10.54 GeV (2.3 fbT 1 ). 



The CLEO III detector includes a high-resolution charged particle tracking system 15] . 



a Csl electromagnetic calorimeter 



161 ]. and a Ring Imaging Cherenkov (RICH) hadron iden- 
tification system [17|]. The Csl calorimeter measures the photon energies with a resolution 
of 2.2% at E = 1 GeV and 5% at E=100 MeV. The tracking system also provides charged 
particle discrimination, through the measurement of the specific ionization dE/dx. 

We detect 7/ mesons through the channel 77' — > 7/7t + 7T~, with 77 — > 77. The branching 



fractions for these processes are (44.5±1.4)% and (39.38±0.26)% 19| respectively. We iden- 
tify single photons based on their shower shape. The photon four-vectors are constrained to 
have invariant mass equal to the nominal rj mass. Subsequently, rj candidates are combined 
with two oppositely charged tracks to form an rj'. Loose 7r consistency criteria based on 



5 



dE/dx measurements are applied to the charged tracks. 

The gluonic rj' production at the T(15) is expected to be dominant only at very high q 2 , 
or, equivalently, at high rf scaled energy Z, where Z is defined as 

ry _ £y _ 2E r)l 



where E v r is the 7/ energy and -Ebeam is the beam energy. Enhanced 7/ production at high 
Z would indicate a large rj'g*g coupling. 

For low-energy rf mesons, photons coming from low energy 7r°s are a severe source of 
background. Thus a it veto is applied for Z < 0.5, whereby photon pairs that have an 
invariant mass consistent within 2.5 o with the nominal tt° mass are not included as the 
candidate photons for 77 reconstruction. We consider only 77' with scaled energy Z between 
0.2 and 1 and divide this range into eight equal bins. Fig. [2] shows the 7/ yields in these bins 
for the T(1S) sample. Fig. [3] shows the corresponding distributions from the continuum 
sample taken at a center-of-mass energy of 10.54 GeV. In order to derive the 7/ signal yields, 
we fit the AM V > V distributions [AM,/^ = M (7r + 77) — M(rf)\ in each Z bin with a Gaussian 
function representing the signal, and a polynomial background. The Gaussian is used only 
to define a ±2.5(7 signal interval. Then the 77' yield in this interval is evaluated counting 
events in the signal window, after subtracting the background estimate obtained from the 
fit function. As the AM V / V signal is not described well by a single Gaussian function, this 
procedure minimizes systematic uncertainties associated with the choice of an alternative 
signal shape. 

Information on the gluon coupling of the 7/ can be drawn only from the decay chain 
T(1S) — > ggg* — > rj'X, described by the Feynman diagram in Figure 1. Thus we need to 
subtract both continuum 7/ production and 7/ from the process T(15) — > 7* — > qq. The 
latter component is estimated using 

B(T(15) ^qq) = R- B{T(1S) -> /j + /0 = (8-83 ± 0.25)%, (2) 

where R is the ratio between the hadronic cross section 7* — > qq and the di-muon cross 



section 7* — > fi + fi at an energy close to 9.46 GeV. We use R = 3.56 ± 0.07 18] and 
B(T{1S) -> (jL+fj,-) = (2.48 ± 0.05)% [19]. The yield N(T(13) -> ggg*) is estimated with 
the relationship 

N(T{1S) ggg*) = iV had - iV( 7 * -> qq) - N(T(1S) -> gg), (3) 
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FIG. 2: The spectra of the difference of the t]tt + tt~ and r] invariant masses in different Z ranges 
reconstructed from T(IS') data, fit with a single Gaussian function for the signal and a first-order 
polynomial for the background. 

where iVh a d is the number of hadronic events in our sample, and N(j* — * qq) is the number 
of continuum events derived from the 10.54 GeV continuum data set, corrected for the lumi- 
nosity difference between resonance and continuum data, and the center-of-mass dependence 
of the cross section for the process 7* — * qq. 

The two dominant components of the rj' spectrum have different topologies: T(IS') — > 
ggg* produces a spherical event topology, whereas qq processes are more jet-like. This 
difference affects the corresponding reconstruction efficiencies. The •ygg cross section is 
only about 3% of the ggg* cross section; thus we make no attempt to subtract the former 
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FIG. 3: The spectra of the difference of the rj7r + Tr~ and r/ invariant masses in different Z ranges 
reconstructed from continuum data taken at a center-of-mass energy of 10.54 GeV, fit with a single 
Gaussian function for the signal and a first order polynomial for the background. 



'or the two event topologies 



201 ] Monte Carlo samples, 



component from the latter. Fig. H] shows the efficiencies obtained 
with CLEO III Monte Carlo studies. We use GEANT-based 
including T(1S") and continuum samples. In order to use the continuum sample taken 
at 10.54 GeV center-of-mass energy for background subtraction, we develop a "mapping 
function" to correct for the difference in phase space and Z range spanned in the two 



samples. The procedure is described in detail in Ref. 



131 ] . In brief, we use the rf energy 



distribution functions for the Monte Carlo continuum samples at center-of-mass energies 



S 



equal to 9.46 and 10.54 GeV and obtain a relationship between the measured Z at higher 
center-of-mass energy Z10.54: 



Z 9A6 = -0.215 x 1(T 2 + 1.2238 Z 10 . 54 - 0.6879 (Z 10 . 54 ) 2 + 0.8277 (Z 10 . 54 ) 3 - 0.3606 (Z 10 . 5 4) 4 • 

(4) 

We derive the 7/ spectrum from the decay T(IS') — ■> ggg* — > t/X from the efficiency 
corrected 77' yields in each Z bin, subtracting the contributions from continuum and T(IS') — > 
gg as shown in Equation^ In this case, the r/ spectrum from the process T(15) —>■ qq —>■ rj'X 
is corrected for the distortion introduced by initial state radiation (ISR) on the continuum 
7/ energy spectrum used to account for this effect. The partial yield N(T(1S) — > qq)(Z) is 
estimated with the relationship 



N(T(1S) - gg)(Z) = X( 7 * - gg - rfX)(Z) x i? ISR x ^H£i 

= X( 7 * - gg _> r/X)(Z) x i? ISR x ^h^I ? (5) 



where -Risr accounts for the difference between the T(1S) — > gg — > r/'X and the 7* — > 
qq — > rj'X spectra due to initial state radiation (ISR) effects, estimated using Monte Carlo 
continuum samples with and without ISR simulation, and a(T(lS) — > qq)/a(e + e~ — > gg) 
accounts for the relative cross section of these two processes. The correction factor -Risr 
differs from 1 by a few percent at low Z and as much as 25% at the end point of the 
7]' energy. The cross sections used are aCT(lS) — > At + /i~) = 0.502 ± 0.010 nb and 



<j(e + e — > fi + fi ) = 1.372 ± 0.014 nb 21(. Fig. [5] shows the measured differential cross 



sections da^/dZ for the processes T(1S) — > ggg*, T(1S) — > gg, and, 7* — > qq. 

Theoretical predictions give the energy distribution function dn/dZ = [1/N (¥(]£) — > 
ggg*)] x dN(T(lS) — ► ggg* — > rj'X)/dZ; we obtain the corresponding experimental quantity 
by dividing by the total number of N(T(1S) — > ggg*), estimated by applying Equation [3] 
without any Z restriction. Figure 0a) shows the T(15) — ► ggg* — > r/'X energy distribution 
function, whereas Fig. [6jb) and c) show the corresponding distributions for T(1S) — > gg — > 
t/'X, normalized with respect of the total number of T(1S) — > gg and T(IS') — * t/'X, 
normalized with respect to the total number of T(15*). 

The inclusive 7/ production at the T(IS') is expected to be dominated by the transition 
T(15) — ► ggg* — > rj'X only at high 77' energy. The energy at which this occurs cannot 
be predicted from first principles: an empirical criterion is the x 2 of the theory fit to the 
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FIG. 4: The rj' reconstruction efficiencies as function of Z for different MC samples with no ir° 
veto, and with ir° veto in photon selection. The 7r° veto was applied in this analysis for Z < 0.5. 
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FIG. 5: The measured differential cross sections daw/dZ for a) 7* 
T(IS) -^qq^ rfX (white) and c) T(1S) -» ggg* -> r/X (black). 



— > r/X (shaded), b) 



data. For example, a numerical analysis of the CLEO II data 14j obtained a x 2 of 2.4 
for three degrees of freedom, using the 3 experimental points at Z > 0.7, and ~ 24 for 4 
degrees of freedom using the 4 points at Z > 0.6. This observation led Ali and Parkhomenko 
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FIG. 6: The energy distribution function dn/dZ as defined in the text for a) T(15) — > ggg* — > rj'X, 
b) T(15) -» qq -> rfX, and c) T(15) -> rfX. 
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to conclude that the Z region likely to be dominated by T(IS') — > ggg* — > rj'X starts at 
Z = 0.7. Thus we quote global branching fractions for T(IS') — > rj'X and the corresponding 
results for Z > 0.7 separately. 

Tabled] summarizes the dominant components of the systematic uncertainties. The overall 
relative errors on the rj' branching fractions are ±8.1% for qq — > rj'X, ±9.1% for ggg* — > rj'X 
for Z > 0.7 and ±7.2% for all other branching fractions. 



Sources 


ggg* Sample (Z > 0.7) 


qq Sample 


All others 


Reconstruction efficiency of ^ 


2.0 


2.0 


2.0 


Reconstruction efficiency of r\ 


5.0 


5.0 


5.0 


Number of rf from fit 


1.0 


1.0 


1.0 


Total number of T(15) 


1.0 


1.0 


1.0 


B(rf — > ir + ir~r]) 


3.4 


3.4 


3.4 


B(T(15) - qq) 




3.0 




Ratio of integrated luminosity 


1.9 


1.0 






0.7 


1.3 




7T° veto 




1.7 


0.4 


Z mapping 


6.0 


3.0 


3.0 


Total 


9.1 


8.1 


7.2 



TABLE I: The components of the systematic errors (%) affecting the branching fractions reported 
in this paper. 



Thus we obtain 
n(T(15) - (ggg* 



rfX) = 



n(T(\S) - (qq) -> rfX) = 



n(T(lS) -> t/X) 



iV(T(lg) -> ggg* -> r/X) 
X(T(1S) - g<?<?*) 
X(T(1S) -> gg -> r/X) 



X(T(l£) - gg) 
_ N(T(1S) -> r/X) 



(3.2 ±0.2 ±0.2)%, 
(3.8 ±0.2 ±0.3)%, 
(3.0 ±0.2 ±0.2)%. (6) 



X(T(LS)) 

The T(IS') — > r/X branching fractions at high momentum (Z > 0.7) are measured to be 

n(T(lS) - (<?g<7*) - t/X) z >o.7 = (3.7 ± 0.5 ± 0.3) x 10" 4 , 
n(T(lS) -> (gg) -> r/X) z>0 . 7 = (22.5 ± 1.2 ± 1.8) x 10~ 4 , 

n(T(lS) -> r/X) z>0 . 7 = (5.1 ± 0.4 ± 0.4) x 10~ 4 . (7) 
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III. COMPARISON WITH THEORY AND CONCLUSIONS 



A. Kagan [6| used the ratio Rz>o.7, defined as 

nth 



R 



Z>0.7 



n 



exp 



(8) 



Z>0.7 



to obtain a first rough discrimination between form factors having drastically different q 2 
dependence. At the time that Ref. ||| was published, the comparison was based on 90% C.L. 
upper limits on the data. This test repeated with our present data give values of Rz>o.7 ~ 74 
for a representative slowly falling form factor Rz>o.7 ~ 25 for the intermediate form factor 
studied by Ref. [5J], and Rz>o.7 ~ 2 for the perturbative QCD inspired shape. Thus the last 
shape is the closest to the data, but it does not match them very well. 

Several perturbative QCD calculations of this process exist, and are based on different 



choices of the form factor H(q 2 ): Kagan and Petrov [6fl assume H(q 2 ) m const m 1.7 GeV 1 ; 
Ali and Parkhomenko relate H(q 2 ) to the expansion of the two light-cone distribution am- 

n 

plitudes (LCDA) describing the quark and gluon components of the n' wave function |14| |. 
Figure [7] shows the measured dn/dZ distribution, compared with three representative choices 
for H(q 2 ): H(q 2 ) = Hq = 1.7 GeV -1 , H as , based on the asymptotic form of the n' meson 
LCDAs, and H(q 2 ) corresponding to the Ali and Parkhomenko Ijj] formalism, with the 
parameters extracted from the previous CLEO II data and the constraints from the rf — 7 



transitions 



111 ] . Note that in the range where T(15) — ► ggg* is the relevant source of rf 



most of the discrepancy between theory and data occurs in the Z = 0.7 bin. In fact, the 
X 2 for the fit of the new data with this theoretical parametrization is 27 for 3 degrees of 
freedom. This may imply that higher order terms in the QCD expansion need to be taken 
into account, or that the T(IS') — > ggg* is not the dominant source of rj', at least at a scaled 
energy as high as Z = 0.7. 

In conclusion we have measured the energy spectra of the rf meson in the decay T(15) — > 
r)'X. Our results are not very well described by existing models based on strong gluonic 
coupling of the rf . Thus the observed B — > r)'X inclusive branching fraction is unlikely to 
be explained by an enhanced g*gr]' form factor, and an explanation outside the realm of 
the Standard Model or an improved understanding of non-perturbative QCD effects may be 
needed to account for this large rate. 
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FIG. 7: Energy spectrum of the r/-meson in the decay T(1S) — > rj'X (open squares correspond to 
the data presented in this paper, filled circles are the previously reported CLEO II data [13]): a) 
measured spectra (log scale); b) expanded view of the Z > 0.7 region to show the comparison with 
the theoretical predictions more clearly (linear scale). The dashed curve corresponds to a constant 
value of the function H(p 2 ) = Hq ~ 1.7 GeV -1 , and the dash-dotted curve (i? as (p 2 )) corresponds 



to the asymptotic form of the r/-meson LCDA 



14|] (i.e., B^ q) = and B { 2 9) = 0). The spectrum 



with the Gegenbauer coefficients [14] in the combined best-fit range of these parameters is shown 

in the shaded region. The solid curve corresponds to the best-fit values of the parameters in the 

form factors from Ref. [ijj] from the analysis of the T(IS') — > rj'X CLEO II data alone. 
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